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% : Abstract 

\ In this paper we expand our previous investigation of a quantum particle 

o, 

subject to the action of a random potential plus a fixed harmonic potential 
at a finite temperature T. In the classical limit the system reduces to a well- 
known "toy" model for an interface in a random medium. It also applies to 

Sh , a single quantum particle like an an electron subject to random interactions, 

O 
O 



where the harmonic potential can be tuned to mimic the effect of a finite box. 
Using the variational approximation, or alternatively, the limit of large spatial 

H : 

dimensions, together with the use the replica method, and are able to solve the 

model and obtain its phase diagram in the T — (h 2 /m) plane, where m is the 

particle's mass. The phase diagram is similar to that of a quantum spin-glass 

in a transverse field, where the variable h 2 /m plays the role of the transverse 

field. The glassy phase is characterized by replica-symmetry-breaking. The 

quantum transition at zero temperature is also discussed. 
05.30.-d,05.40.+j,75.10.Nr 
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I. INTRODUCTION 



The effects of quantum ffuctuations on phase transitions is a topic of significant current 
research. On the other hand transition into a glassy phase in disordered systems is a topic 
which is far from trivial and requires special analytical and numerical techniques. In par- 
ticular the behavior of interfaces in disordered systems has been the subject of many recent 
papers [1-12]. It is of interest to investigate the combined effect of disorder and quantum 
fluctuations that play a significant role at low temperatures. 

In a recent paper [JHJ , we have focused on the problem of a quantum particle in a random 
potential with power law correlations plus a fixed harmonic restoring force. This problem 
has previously been studied classically |li|8|,|9 ,[11,121, by using the variational approximation 
fi]9[]Tllj (or alternatively the large iV-limit IflOj, 1 2 1 , where N is the number of dimensions). 
A transition into a glassy phase has been identified. The signature of such a transition is 
the appearance of solutions with replica-symmetry-breaking (RSB) ||14|| . RSB is typically 
associated with a complicated free-energy landscape characterized by many local minima 
which are separated by high barriers. A sharp transition for a single particle must be an 
artifact of the variational approximation ( or alternatively the large- N limit). In reality this 
must be a crossover effect over a range of temperatures. 

The "toy" model of a particle in a random potential in one dimension is the simplest 
example of an interface problem where the location of the particle stands for the location of 
an interface between to phases (like up an down spins) and the random potential stands for 
quenched impurities which can pin the interface. Understanding of this problem is the first 
step of investigating higher dimensional manifolds. But the quantum version of the model 



also apply directly to a particle like an electron in a dirty metal [17], where the strength of 
the harmonic potential can be tuned to mimic a finite box which the electron is confined to, 
like the situation which occurs in certain mesoscopic systems. 

The question we have had in mind was what were the effects of quantum fluctuations on 
this glassy transition. There are some similarities between this problem and the quantum 
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spin-glass problem in a transverse field In that problem, tunneling effects produced 

by the transverse field eventually lead, for strong enough field to the destruction of the 
spin-glass phase. But this phase when it exists, is still characterized by replica-symmetry- 
breaking (RSB), for the infinite range model. With this analogy in mind, one should notice 
though, that the transition which has been found for a particle in a random potential P,P|JTT| 
is of the Almeida-Thouless type in the sense that it is associated with the appearance 
of RSB but not with an order-disorder (spin-glass like) transition. 

In this paper we expand our previous investigation of the quantum particle, and also 
correct a typographical error that unfortunately occurred in equation (7) of |13|| , propagated 
into eq. (8) and affected some of the numerical results. We find as before that the glassy 
phase, characterized by RSB persists in the presence of quantum fluctuations. However, 
we find that for small enough particle mass, (m and H enter only in the combination k = 
(m/h 2 )) the glassy phase ceases to exists. Thus the parameter 1/k plays a similar role to 
the transverse field in the Ising spin-glass. A schematic phase diagram is depicted in Fig. 1. 



II. BASIC DEFINITIONS 

The density matrix for a particle at finite temperature T = l/ks/3, subject to an har- 
monic potential and a random potential V, is given by the functional integral JTJ] : 

r x([/)=x' f i 



rx(U )=x- I r 

p(x,x',C/) = / [dx]exp{--/ 
ix(o)=x I n Jo 



m±(uf + tMuf + } 



du}, (2.1) 
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where x is a A-dimensional vector (N is the number of spatial dimensions), and U = 13%. 
The variable u has dimensions of time and is often referred to as the Trotter dimension. For 
the diagonal elements of p we observe that the trajectory x(w) forms a closed path. In this 
paper we are concerned with a random quenched potential V(x) and quantities of interest 
are the averaged free energy and the mean-square-displacement given by 

- 13(F) a = (\n J p(x,x) dx), (2.2) 



= 7p(xx)x^x 
\ J P(x, x) dx I R 



where, Qr stands for an average over the random realizations of the potential. The difficulty 
in carrying out the quenched averages stems from the fact that in equation ( |2.2| ) the average 
is taken after taking the logarithm, i.e. the averaged free energy is not the logarithm of the 
averaged partition function. Similarly in eq. ( |2.3| ) the average of the quotient is not equal 
to the quotient of the averages. We take the potential V(x) to gaussian distributed, which 
means that the probability for a given realization of the potential is 

P(V(x)) = Cexp(- J dxdx'V(x)A(x - x')V(x')), (2.4) 

with some known function A(x — x'). It is thus sufficient to know only the first two moments 
of the distribution vis. 

(v(x)) H = o, <ywvw) B = -Nf (^4r^) ■ (2 ' 5) 

where the functions / and A are related to each other. The function / describing the 
correlations of the random potential is taken to decay as a power at large distances: 

f(y) = (go + y) ^ (2.6) 

The index 7 describes the behavior of the correlations of the disorder at large distances. In 
this paper we consider only two cases, that of 7 = 3/2 which we call the case of short-ranged 
correlations, and that of 7 = 1/2 which we call the case of long-ranged correlations [|7||8|,pT| . 
The parameter do plays the role of a short-distance regulator for /. We chose these values 
of 7 in order to make contact with known results in the classical case. In the classical case 
it has been shown (for N = 1), that faster falling correlations than those characterized by 
7 = 3/2, even correlations falling exponentially fast, are equivalent within the variational 
approximation to the case of 7 = 3/2 |7yTT|j. This fact also holds in the quantum case, as is 
demonstrated in sec. Ill below. The long-range case is of interest because of its connection 
with the directed-polymer problem and with the random-field Ising model P,|7|j9|JTl| . 

We apply the replica trick in order to carry out the quenched average over the random 
realizations. We consider n-copies of the system, and obtain for the averaged density matrix: 
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rX a (U)=X a " 

p(xi ■ ■ -x n ,xi ■ ■ -x n , U) = / l[[dx a }exp{-H n /h} , (2.7) 

Jx„(0)=x a 0=1 
1 /•[/ f , 



Note that this equation although somewhat similar, is quite different from the equation for 
the n-body hamiltonian corresponding to the directed polymer problem in 1 + N dimensions 

0. The difference is first the fact that the integral over the "time" variable u is on a finite 
interval, and second and more significantly, that the n-body potential is non-local in "time" , 

1. e. it involves a double integral over both u and u' . 

III. THE LARGE-iV LIMIT AND THE VARIATIONAL APPROXIMATION 

In the large- N limit, we introduce collective variable fields r ab (u, u'): 

r ab (u, u) = jj*a{u) ■ xt(u') (3.1) 



and Lagrange multiplier fields s ab (u,u') to implement (|3.1|) . We find fllO 

rX a (U)=X a ™ r 

p(xi • ■ ■ x n , x x • • ■ x n , U) = / [j [dxa] / j| [dr ab (u, u')\ [ds ab (u, u)) 

Jx a (0)=X a a=1 J a < 6 

exp{-— 2/ du du' 2_,( Nr ab (u,u') - x a (u) • x 6 (w')) s ab (u,u') 

1 f u r 1 
~2h Jo du ^ i mx ^ M ) + W 

-72// dudu'^Nf (r aa (u,u) + r bb (u',u') - 2r ab (u,u'))\ . (3.2) 

ih J Jo ^ J 
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From here it follows that the free energy is given by 
1 r u r u 



j ru ru 

,(3{F) R /N = —2 du du'y2(r ab (u,u')s ab (u,u') + 
2h Jo Jo ^ 

/OO /* 
Y[dx a / f|[(ia; a ]e J , (3.3) 

-00 „ Jx„ (0)=x„ (U)=x„ 



where we defined 



1 r u / \ 1 r r u 

J= ~2hJo du ^ ( mi a( M ) +^a( M )J + 7^2 J J o du du'J2 s ab(u,u')x a (u)x b (u'). (3.4) 



The limit n — > is to be taken. In equation (|3.3| ) the variables x a (u) are scalers since a factor 



of N has been extracted. In the large N limit the collective variables r a b(u, u') and s a b(u, u') 
are determined by the stationarity of the free-energy. We will be looking for solutions of 
the saddle-point equations obeying translational invariance in the "time" direction, i.e. such 
that the order-parameters depend only on the difference ( = u — v! . They should also be 
periodic functions of this variable with period U, and also symmetric under ( — + —( (even 
functions). 

In the case of a classical particle there is another approach to the problem called the 
variational approximation which yields essentially the same results as the large-iV limit 
aside from a simple renormalization of the correlation function of the disorder. For the 
quantum potential we were able to show that the same results still hold. Let us introduce 
the variational hamiltonian: 

1 f U 



I r u r 
h n = - / du pxM + /^( 

^ Jo a 

1 f U f U 

— -zt / du du'y^SabOu - u')x a (u) ■ x 6 (u')- (3-5) 
2n Jo Jo ~f 



ab 

Then the variational free-energy is given by: 



n/3 (F) R /N = {H n - h n ) h Jh - In / [dx]e~ h "/ h . (3.6) 



In Appendix A we show that the variational free-energy coincides with that given in 
the large- N limit (see eq.( |4.5| ) below), with the following renormalization of the function / 
defined in eq. ( |2.6| ): 

The effective correlation / is similar to / in the sense that for large arguments: 
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f(v) ~ xf^V 1 ' 1 ' (3-8) 



with 



N 

7 = 7 if 7<y + l 

N „ N , , 

7 = y + l if 7>y + l (3.9) 



r(l-7 + JV/2) /iV\- 1+7 
9 = 9 wL/on • ( 3 - 10 ) 



and 



r(iv/2) V 2 

Since we will be mainly concerned with the case where 7 = 7, the the values for g quoted 
in the numerical work below should actually pertain to the renormalized coupling g. In the 
limit N — > 00, / approaches /. 

IV. THE STATION ARITY CONDITIONS AND THEIR EXACT SOLUTION 

ASSUMING REPLICA SYMMETRY 

In this section we will derive the stationarity equations for the free energy and solve 
them in the replica-symmetric (RS) case. The case of RSB will be discussed in the fol- 
lowing section. We proceed by transforming to frequency space by using a Fourier series 
representation of the order parameters: 

00 

r a b(C) = -q J2 exp(-iwj C) r ab (ui) (4.1) 

r 1 l = — oo 

00 

Sab(C) = -3 exp(-iwi C) s ab (uji), (4.2) 

P l=-oo 

where 

2n 

uJi = —l, Z = 0,±l,±2,... . (4.3) 

The last equation follows from the periodicity of the order parameters. From the fact that 
^ab(C) (and also s) are even functions it follows that r a b(u>) = f a b(—uj), and similarly for s. 
To proceed we note that for any periodic function g(u — u') with period U it follows that 
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rU rU nil ru rU rO 

du du'g(u -u')= du d(g(() + du d(g(() = 

Jo Jo Jo JO JO J-U+u 

[ U du I" d(g(() + f U du [ U d(g(( - U) = [ U du f° d(g(() = U f dQg(Q. (4.4) 
Jo Jo Jo Ju Jo Jo Jo 



In terms of the new variables introduced in eq. Q4.2|) , the free energy becomes: 

n(3(F) R /N = ^EE^W-^Etr lnG M 

Z U) ab Z w 

+4 E f d Cf(lj2 (f aa (u') + f w (u/) - 2e~^ f^o/))] + const., (4.5) 

ab ® \P uj' ) 

where we defined 

Gab(uj) = ([(m^ 2 + - 5(0;)]-%. (4.6) 



The constant in eq. (|4.5|) can most easily be determined from the known free energy for a 
free particle (/ = and u = 0). 

From eq. ( f4.5|) it follows that stationarity equations are given by: 

f ab {uj) = G ab (u), (4.7) 
~s ab (uj) = \ £ d( e^f (± £ + f»(«/) ~ 2e-^f a6 (o/))) a^b, (4.8) 

5 aa (o;) + £ 5 a6 (0) + lj U dC(l- e^)f{- £ f aa (o/) (1 - e~^)) = 0. (4.9) 

bjta " P ui' 

In eq. Q4.8|) , /' stands for the first derivative of the function /, characterizing the correlations 
of the random potential. 

We proceed to look for a RS solution to these equations by taking all the elements 
r ab (u) ee r(u) with a ^ b to be equal to each other (for the same frequency) and similarly 
for s ab (uj) ee s(lj). The diagonal elements are also taken equal to each other, and are denoted 



by Tdiyj) and 5<j(o;). In the limit n — > we find from equations (|4.9| ), ( |4.8| ) by inverting an 
n x n matrix: 

1 s{uj) 

fdiuJ) = 7, Z , \ ZTT + 7 o ~ / \ ~i wo ; (4-10) 

(muj z + a - s d {uj) + s{uj)) z 
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These expressions have to be substituted into eqs. ( |4.8| ) and ( |4.9|) to obtain the self-consistent 
equations for s and s d : 

m = 2 -[ u dC r (- y ( - + g>0 (1 = e ~^ C) )) (4 12) 

y ' hJo * J \{3j? W 2 + /J - s d (u') + s(u') (muj 2 + Li-~s d {u}) + ~s{uj)) 2> ) ' K ' 

2 r u 



s d {u>) = 5(0) - - / d( (1 - e**) 



fV^. V (1 — p~ iu> '<>) - i 1 (4 TV) 

J mu' 2 + li - S d {u/) + s(oj>) (mu' 2 + /i - s d {u/) + s(u')) 2> - 1 ' ' 

We will now seek a solution with time-independent off-diagonal elements. This is in analogy 



with quantum spin-glasses |L5],[l6|, the rational being that the off-diagonal elements of r 
and s are the glass order parameters and hence are constants independent of the Trotter 
time. We find that such a solution satisfies the saddle point equations exactly. The diagonal 
elements are still time-dependent. In frequency space we look for solutions of the form: 

s{u) = s S u , . (4.14) 



This ansatz satisfies eqs. (|4.12| ), and we obtain: 



''-^(/Hs ^+zUm ) (415) 



2 fU 
s d (u) = s - - 
n 

Two physical quantities of interest are 



in i oo 

((* 2 ))r/N = - X>«(°) = r„(0) = - £ ?i(f*k), (4-17) 

71 a=l P fc=-oo 

in in 

((x 2 ) - (x) 2 ) R /N = - £r aa (0) - -, 77 E r «*(°) = ^(0) - r(0). (4.18) 



, n(n — 1) 

a=l v / afo 



The very last equality holds if one assumes replica symmetry. 

« x2 > ViV = ±- + \ £ * z-r-. + -f~ 2 (4.19) 



and 



((x 2 )-(x) 2 ) R /iV=-^ + i V 1 . . . (4.20) 



In all the equations derived above, m and % enter only through the combination 



771 

K=^i . (4.21) 

h 2 



The classical limit is given by h — > and consequently « — > oo. Practically, this limit can 
be achieved by taking the particle's mass to be very large. On the other hand, the quantum 
limit is obtained for very small particle mass. In the classical limit only the zero frequency 
component of the observables survives and we obtain e. g. 

«^«^+M£)- (422) 

which agrees with refs. [jsyil"]. 

From equations (|4.17|) and ( |4.16 ) it follows that the expression for ((x 2 ))r in the quantum 



case, is the same as in the classical case but with a renormalized temperature 1/(3r given 
by: 

1 1 1 1 1 1, ,„ / 

-T = —r + ~R E 7TZ r T~^ = ~R + o 6 " & K ' ^' 9 > ( 4 - 23 ) 



where we defined: 



bo(/3, K,^g) = W 1 (4.24) 



The minimum of the mean-square-displacement in eq. (|4.22 ) is attained for 



Tf L = \jf c L = ^(2 7 g)WrfT3 - ia ^ (4.25) 



In the quantum case, it follows from eq.( 4.23 ), that the temperature for which the minimum 
is attained, is given by a solution of the equation 

T c + ^&o(l/T c , k, /i, </) = Tf L (4.26) 

We will see in the next section that for the case of long-ranged correlations the temperature 
for which the minimum of the mean-square-displacement is attained is indeed the transition 
temperature into the glassy state characterized by RSB. For the case of short ranged cor- 
relation of the disorder, even in the classical case the transition into a replica-broken phase 
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occurs only approximately at the temperature given by eq. ([4.25|) if \x is not too small [IT . 



Otherwise the transition T* into an RSB phase occurs at a higher temperature than given 
by ( [4.251 ). The value of the mean-square-displacement at the minimum is given by 



((x 2 )) R /N = (1 + 7 )(2 7 )-T^( ? ^ yU -T^ _ i ao . (4.27) 



We now turn to the numerical solution of equations ( [4.161) . Some technical details of 
the numerical procedure are given in Appendix B. First we consider the case of long-ranged 
correlations of the disorder with 7 = 1/2. In Fig. 2 we display the results for the mean- 
square-displacement for various values of k = m/h 2 . We have chosen g = 2v^2, c*o — 0.01 
and fj, = 1 in order to compare with results of the classical case |J. For these values of the 
parameters, T c (oo) = 0.995. For large k there is good agreement with known results for 
the classical particle ||. As k decreases and quantum effects become more significant, the 
minimum of the curves moves to the left which shows that bo is increasing and hence T c is 
decreasing (see eq. (|4.26|) . For k ~ 0.17 the minimum of the curve is at about T = 0, and 



for lower values of k no minimum occurs at a positive temperature. From this data we can 
reconstruct an approximate phase diagram in the T — 1/k plane, see Fig. 3. 

We now turn to the case of short-ranged correlations of the disorder, i.e. 7 = 3/2 
(recall eq. ( p.6|) ). We have chosen ao = 1, g = 1/V2n and /i = 0.53 in order to compare 



with the results of the classical case |TTJ. In Fig. 4 we plot the mean-square-displacement 



given by eq. (|4 . 1 9|) for various values of k as a function of t = T/T c (oo), where T c (k) is the 
temperature for which the "transition" into the glassy state occurs. Thus t=l is the rescaled 
temperature for which the transition occurs in the classical case (k = oo). For the values of 
the parameters we have chosen, T c (oo) = 0.2082. We see that the situation is similar to the 
long-ranged case. The value of k for which the minimum of the curve hits zero is about 2. 
The reconstructed phase diagram is given in Fig. 5. 

In Appendix III we review, correct and present more details of the real-space approach 
used in our previous investigation [13]. The formalism is set up for first-step RSB. (The 
need for RSB at low temperature will be discussed in the next section). We have checked 
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numerically that for replica symmetry and long-ranged correlations of the disorder the plots 
of the mean-square displacement vs. temperature coincide with those obtained from the 
momentum space calculations as reported above. 

We finish this section with a discussion of the zero temperature limit. In this limit the 
frequency becomes a continuous variable, and the equations for § and Sd(u) become: 

§/P = 2f'(- ' 



IT J-oo KUj' 2 + fM — Sd(u!'/h) 

/oo 
dC(l - e^ c ) x 
-oo 

1 r°° 1 — p- iuJ '^ 1 r°° du' 



7T J-oo KUJ' 2 + fl — Sd(00'/h) 7T J-oo KUj' 2 + fi — Sd{oj' /Tl) ' 

s d (0) = 5- (4.28) 

where we have rescaled the variables ( and uo by 1/h and % respectively. In terms of the 
solution to the last equation, the expression for the mean square displacement become: 

((x 2 )WiV = ~&o03 = oo, k, fi, g) + 2/'(& (/3 = oo, k, h, g)). (4.29) 

with 

1 r°° du 1 

bo(P = oo,K,fi,g) = - — — _ (4.30) 

7T J-oo KUJ + H — Sd{Uj'/n) 

Eq. ( |4.29|) is similar to the classical expression ((4.19[) , with the temperature variable T 
replaced by (l/2)yu6 (/9 = oo, k, /x, g). From eq. (|4.26| ) it follows that the quantum transition 



at T = (for the case of long-ranged correlations) occurs when k = k c , where 

1 rpCL 

-b (P = oo,K C7 n,g) = -2—. (4.31) 

2 n 



V. REPLICA-SYMMETRY-BREAKING SOLUTION. 

In this section we proceed to explore the possibility of RSB. The need to break replica 
symmetry has been demonstrated in the classical case below a certain transition temperature 
T c . The physical origin of RSB is the existence of many local minima of the free energy, 
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separated by barriers. A sharp transition temperature for a single particle exists only in 
the large-iV limit or within the framework of the variational (gaussian) approximation. 
Classically the pattern of RSB depends on the range of correlations of the random potential, 
for short range correlations a one-step RSB has been found sufficient, whereas for long- 
ranged correlations a continuous RSB has been found necessary P,|TT||. In the quantum case 



we will show below, that RSB also occurs if the strength of the quantum fluctuations is not 
too large. Support for this assertion also comes from the fact that for the RS solution a plot 
of {{x 2 ))r vs. temperature (Fig. 2) is not monotonic, but increases as T — > for k not too 
small. This indicates that the RS solution is inadequate at low temperatures in this case. 
Practically one can compare the free energies associated with the RS and RSB solutions (if 
the latter exists) and verify that the free energy of the RSB solution is higher ( an artifact 
of the n — > limit). This indicates that RS has to be broken. This is indeed the situation 
in the classical case and it carries over to the quantum case as explained below. 
In constructing a RSB solution we associate 

s aa (uj) = s d (u), (5.1) 
s ab (uj, z) <-> s(z) 8 ufi , (5.2) 

where the Parisi parameter < z < 1 labels the 'distance' between the replica indices ab. 
Again, we have used the static ansatz for the non-diagonal matrix elements, which will turn 
out to be a consistent solution. One should not confuse the frequency dependence of the 
diagonal element §d(u) with the z-parameter dependence of s(z). A similar parametrization 
applies to r a b(oo). To write down the expressions for f, analogous to eqs. (|4.11|) , we need the 
expression for the inverse of a hierarchical matrix (see Appendix II of ||). We then find: 

~ / \ 1 A , s{z = 0)(L,o . f 1 dz [5](z) \ 

muj 2 + ix - s d (u) + (L,oJo dz s ( z ) V ^ Jo z f x +[ s \i z )J 

~, \ ~, c 1 ( [*K Z ) s(z = 0) f z dz \sMz) \ , r A . 
r(u,z) = r{z)5 Ut o = 5 Ui0 -[ L n ; + ^ + / — l i } J , 5.4 

fJL \Z{fJL + [S\(Z)) LI JO Z 2 IX+[S\(Z)J 

where 
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s](z) = z s(z) — dz s(z). (5.5) 







Using these formulas, the stationarity equations (see eqs. (|4.8|) and (|4.9|)) become: 

S(Z) = Wf {p [z{fi+[8\{z)) ~ L~z* »+[s)(z)) + P 5 m ^ 2 + H ~ ' (5 ' 6) 



f dz ~s{z) - I r d( (1 - e **) f [W (1- e-^) 
Jo hJo \Pj^o 



S d (u>) = I dz ~s(z) -~-\ dC (1 - e«") f | - R £ (1 - e-"<) 1 ] . (5.7) 

If we compare equation ( |5.6|) with the corresponding equation in the classical case, we 
see that it is the same , apart from a "renormalization" or shift of the constant ao appearing 
in the definition of the function / (see eq. (|2.6|)): 

ao -> clr{/3, m, /i, g) = a + b (/3, k, /i, g). (5.8) 

The crucial point is that b is independent of z. We will discuss first the case of long-ranged 
correlations of the disorder. One can repeat the steps carried out for the classical case fl8l,|TT 
and we find the following solution to eq. (|5.6|): 



slz) = < 



\Az\ < z < Z\ 



3 
2 



Az z\ < z < z 2 

-^Az\ 



%Az\ z 2 < z < 1 



with 



A = (2/3)V/3 3 (5.9) 
3 



Zl = 2 9 

and Z2 is the solution of the equation 



-2/3, ,1 



/i 1 / 3 /?- 1 (5.10) 



±PAa R 4 + Z2-^ = 0. (5.11) 

We observe that the only difference from the classical case is the appearance of the 
renormalized function a# instead of ao in eq. ( |5.11| ). This solution becomes replica symmetric 
for Z\ — z%, which occurs when 
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T c = 1//3 C 



I^V 1/S 



-a R {(3 c ,m,fi,g) 



(5.12) 



This constitutes an equation for the transition temperature into the glassy phase. Notice 



that for 7 = 1/2 this is the same as equation ( 4.26 ) in the previous section, which was 
obtained by minimizing the mean-square displacement in the replica symmetric case. Also, 
it follows from eqs. ( |4.17| ) and ( |5.3| ) that in the presence of RSB: 
2w , AT 1 s(z = 0) 1 f 1 dz [s](z) 



_/ ' 1 1 

Pn Jo z 2 fi+ [s] (z) p ^ mu 2 s d (u) 



(5.13) 



:9 ' V 



(5.14) 



Using the RSB solution obtained above, we find: 

° r -'ir'-'-' (I,, - -hn = - . 

2" ^ 2 

which is exactly the same as found for the minimum of the replica symmetric expression 
in the previous section (eq. ( |4.27 ) for 7 = 1/2). Thus, through the entire glassy phase the 
mean-square displacement remains locked at this constant value. 

Let us comment briefly on the zero temperature limit of the RSB solution. In that limit 
the solution given in eq. ( |5.9| ) becomes 

(gn) 2 ^ < z < z x 

(l) 2 (gl3z) 2 z 1 <z<z 2 , 

ga R (oo, m, /i, g)~ 1/2 z 2 < z < 1 



slz 



P 



with 
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- l ' 2 a R {oo,m.^g)- l ' A p-\ 



(5.15) 
(5.16) 



For P — ► 00 this solution appears more replica-symmetric like but it really isn't, since the 
value at z=0 is always different from the value at z 7^ and the contribution to the mean- 
square-displacement remains flat. 



In the short-ranged case equation ( |5.6| ) has a one-step RSB solution, as in the classical 
case, again since the only difference from the classical equation is a renormalization of a®. 
The equations can only be solved numerically |TT| , and a detailed investigation will be carried 
elsewhere. 
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VI. CONCLUDING REMARKS 



In this paper we have investigated the combined effects of quantum fluctuations and 
quenched disorder for the case of a particle subjected to a combination of a fixed harmonic 
potential and a random potential with power law spatial correlations. Acting separately, dis- 
order and quantum fluctuations both increase the mean-square- displacement of the particle. 
From Figs. 2 and 4 it becomes evident that at high temperature adding quantum fluctua- 
tions in addition to disorder increases the mean-square-displacement but at low temperature, 
in the phase characterized by glassy behavior, the mean square displacement remains locked 
probably due to tunneling effects among the different free-energy minima. 

We have seen that the phase diagram of the model is similar to that of the quantum 
spin-glass in a transverse field. There is a transition at zero temperature from a phase 
characterized by RSB effects to a phase where such effects are not present. 

It may be possible to extend the problem of a quantum particle, to higher-dimensional 
manifolds in a disordered medium, and explore the effects of quantum fluctuations. It may 
be also possible to extend the problem to many particle systems where statistics might play 
an important role. We hope our work will stimulate further research in this interesting area. 
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APPENDIX A: 



The major step in calculating the variational free-energy is evaluating the expectation 
value: 

-*w>' )) t , (A1) 

This is done by expanding / in a power series and using the following formula: 



((XoM - XkM)*)^ = (JV - 2)!! ' G °°<°> + Gi *(°) " 2G <*(" " "'O 3 • < A2 > 
where we defined 

iVG a6 ( M - u') = (xa(u) • x^u'))^ • (A3) 

Resumming the series gives rise to f(G aa (0) + GWO) — 2G a b(u — u')) . 

In momentum space conjugate to the Trotter time the variational free-energy becomes: 

n(3(F) R /N = WE G a b ^)s ab (co)-Wtr 

Z u ab A w 

+1- E dC ^ (7? E ( G -(^0 + G 66 (u/) - 2e-- '< G ab (u'))) + const. (A4) 

where s ab (u)/j3 is the Fourier transform of s a b(u — it') and G a b(co) os the Fourier transform 
of G a b(u — u'). This free-energy coincides with the one derived from the large- N limit with 
the replacement /—>•/. 

APPENDIX B: 

In this Appendix we provide some further details on the numerical solution of equations 
Q4.16Q and the evaluation of the mean-square-displacement, eq. ( |4.19|) . Our goal is to 



minimize the error in the truncation of high frequencies. Hence we define the function td(u) 
for uj ^ Oby 

i d {u l ) = s d {u) + i (Bl) 
17 



Where t will be adjusted such that for the highest frequency included, t<f(a> r , 
then use the formula plfl: 



with 



to write 



where 



1 ^ e 



1 cosh(a:(l -2C/U)) 1 
sinh(a) 



-Y — 

P J^O mU}2 + V " S d( UJ ) 
A k - 

+7iY 



p m 

2 V K 



2 1 

+ 



cosh(a R (l - 2C/U)) 



e- l ^t d (uj k ) . 

P t^i ( muJ k + A* + t){™ul + fi + i- t d (uj k ) J 



0. We 



(B2) 



(B3) 



(B4) 



Hr = H + t 

aR ~ 2 V K ' 



(B5) 
(B6) 



A similar expression is used for the case of ( = 0. The equations for td(u>k), t and s were 
solved numerically with k max = 10. We used a fortran routine [|TJJ that finds a root for 
a set of non-linear equations. The solution was then used to calculate the mean-square 
displacement. 



APPENDIX C: 



In this Appendix we give some of the details omitted in ref. JT3J for lack of space, as well 
implement some corrections. Some of the notation is slightly different than in the present 
paper. Starting with eq. ( |373| ) we proceeded in (I) to look for a saddle point under the 
constraint that the off-diagonal elements of the order parameters r and s are independent 
of the Trotter time, whereas the diagonal elements denoted by x( u ~ u ') = r d(u — u') and 
v{u — u') = Sd{u — u') are dependent on u — u' . In the 1-step RSB we have used the notation 



r(z) = r 2 , z < k (CI) 
r(z) = r n , z>k (C2) 

where k has been used to denote the breaking point z c . A similar notation apply for s. 
Substituting these order parameters in the expression for the free-energy we obtain: 

(3{F) R /N = — [k(r 2 s 2 - r n s n ) - r 2 s 2 ] + / / dudu'x(u - u')v(u - u') 
Z zh Jo Jo 

_^fc/(2 X (0) - 2r u ) + ^(k - 1)/(2 X (0) - 2r 2 ) 



+ 



wll dudu> ^ 2x ^ ~ 2x<yu ~ u '^~\S°° Dz ln / Dy ( ZH ^ k + const ' ^ 



2h z 

where Dz = (dz/\/2n) exp(— ^ 2 /2) and similarly Dy, and Z H = / f[dx]e~ H / n . H is the 
one-dimensional, one-particle effective hamiltonian given by: 

if = / du —x 2 {u) + -x 2 (w) - (zy/s^ + yV«2 - sn)a;(«) 
Jo l Z Z 

1 

— - / / dudu' (v(u - v!) - s 2 ) x(u)x{v!) (C4) 
zh Jo Jo 

The different order-parameters are determined self-consistently from the equations which 
extremize {F) R : 

X(u-u') = jD, iD ^f^« Z » , (05) 
r2 _ , „_ JDy(x(u)) 2 H Z k H 



f B7 )_uy {*W)h*h (r7] 
J Dz — Td^zI — ' (C7) 



fDyZjfr 

v[u - u') = [2fc/'(2x(0) - 2r n ) + 2(1 - fc)/'(2 X (0) - 2r 2 ) 

2 /-^ /-^ 
"fT 2 



Jo 



dudu'f'(2 x (0) - 2 X ( M - u'))]U5{u - u') + 2/'(2 X (0) - 2 x (u - u')) (C8) 

Sll = 2/'(2 X (0)-2r n ), (C9) 
s 2 = 2/'(2 X (0)-2r 2 ). (CIO) 

The equation for the breaking point k is given by 

y (r 2S2 - r liail ) - ^/(2 X (0) - 2r u ) + ^/(2 X (0) - 2r 2 ) 
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We have continued by putting the "time" variable on a lattice, with lattice spacing U/M, 
where M is to be taken eventually to infinity. In practice we carried out exact calculations 
up to values of M = 20, and extrapolation to M = oo has been made by finite size scaling 
and the 1/M 2 rule The calculation involves finding the inverse and the determinant of 
the M x M matrix M. appearing in the expression for the discretized effective hamiltonian, 
see eq. (|C4|) : 



Mij = Mtej + ^) - jf-faij - -s 2 ) i< j, i ± 1, (C12) 

where k = m/h 2 . For i — 1 there is an additional term {—kM/0)5m,3- For j < i the matrix 
is given by symmetry. 

Using this notation, the final result for the discretized free-energy is: 

P(F)r/N = ^-[k{r 2 s 2 -r llSll )-r 2 s 2 } + £ X (|i - j\)v{\i - j\) 

-M!/(2 X (0) - 2r u ) + {k ~ 2 1)/32 /(2x(0) - 2r 2 ) + ^ E /( 2 *(°) ~ 2 *G* ~ jl)) 

_^ i n (iM.) + I i n det M + -1 ln(l - kav 2 2 ) - , ° V ] * , (C13) 
2 1/3/2 2fc v 2; 2(1 - jfetruf) V ; 



where we defined 



CT = ^E(- M_ %' v? = /?V, ^I = /? 2 (s 2 - Sll ). (C14) 
j 



The different order parameters have been obtained by extremizing eq. ( |C13| ). 
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Figure captions: 

Fig. 1: Schematic phase diagram of a quantum particle in a random potential plus a 
fixed harmonic potential. 

Fig. 2: Plot of the mean-square-displacement vs. t = T/T c (oo) for the long-ranged 
case (7 = 1/2). Curves from right to left are for k = 1000,2,1,0.5,0.25,0.2,0.18,0.1. The 
horizontal dashed line represent the RSB solution below T c (k). 

Fig. 3: Phase diagram for the long-ranged case for g = 2a/2 and /J, — 1. The dashed line 
is a possible extrapolation to T = 0. 

Fig. 4: Plot of the mean-square-displacement vs. T for the short-ranged case (7 = 3/2). 
Notation as in Fig. 2. 

Fig. 5: Phase diagram for the short-ranged case for g = l/y/2n and fx = 0.53. Curves 
from right to left are for « = 1000, 10, 4, 3, 2.5, 2, 1. 



23 




T 



3.5 



3.4 



3.3 



DC 



CM 
X 



3.2 



3.1 




T /T c (oo) 



